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- - ■ Abstract 

(N 

^^^ , The quark-mass dependence of the r] in the Schwinger model, which - hke the r]' in QCD - 

p^ ' becomes massive through the axial anomaly, is studied on the lattice with Nj = 0,1,2. 

^ ■ Staggered quarks are used, with a rooted determinant for Nf = l. In the chiral limit the 

D , Schwinger mass is reproduced, which suggests that the anomaly is being treated correctly. 

o 

(N, 

1 Introduction 

_ Staggered fermions [I] offer a cost-effective way of regulating QCD with four degenerate species. 

(j_) ! In nature, the four lightest quarks, known as the m, d, s, c flavors, are far from being degenerate; 

;-^, | only the u and d quarks are approximately degenerate in the sense that m,d — m,u <^ Ahad, 

where Ahad = 0(1 GeV) denotes a typical hadronic scale. Since for an integer number of 

degenerate dynamical flavors the functional measure of QCD scales with the iV/-th power of 

O ■ the determinant, it has been proposed |2j to reverse this relationship, and to represent Nf = 2-1-1 

/^ . QCD at finite lattice spacing a (or cut-off a~^) by the EucUdean partition function 
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Z = JDU dei^'\D,,^^^^J detl/^(D,,ag,n^J e"^- (1) 

where the path-integral runs over all gauge backgrounds f/, and S^ denotes the gauge action. 
Thus the square-root of the determinant of a staggered field with the isospin averaged light 
quark mass irtud = (^^M + ^^d)/2 and the fourth-root of the determinant of a field with the 
strange quark mass mg are utilized to define the regulated version of QCD which is used in 
several state-of-the-art studies of phenomenologically relevant quantities (see e.g. [3]). 

In recent years, the setup ([1]) has been criticized [1H6], because there is no field-theoretic 
proof that its continuum limit is really QCD, or put differently, that the lattice theory ([1]) is 
in the correclo universality class. The issue is more involved than ([T]) would suggest, since in 
practice one needs the generating functional Z[fi, rj] rather than the partition function, and the 
manner in which a given staggered field is reduced to a single "taste" (the modern word for a 
single species within a staggered field) in the valence sector differs from the rooting recipe ([1]) 
that is applied in the sea sector of the theory. Accordingly, the question is whether these two 
reduction mechanisms work in concert, to define a valid discretization of QCD. 



-^In addition to the summary talks [TlflO], the interested reader is referred to the Schwinger model condensate 
tests of [TT11I3] , the eigenvalue based arguments of |12|141[T7] , the analysis in rooted staggered chiral perturbation 
theory [T81IT9] , the renormalization-group based arguments of [201121] , and the analysis of the 't Hooft vertex [22] . 
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Table 1: Overview oi the eL = 17.89 simulations (top). The entries in the M^L (form = 0) and 
MtjL columns quote the predictions for Nf =2 f23 |, [23| / . For Nj =1 only the former exists and 
is smaller by a factor \/2. To test for finite volume effects, the (3 = 3.2 runs are repeated at 
eL = 13.42, 22.36 (bottom). Simulations are performed at Nf =0, with reweighting to Nf =1, 2. 



As much of the criticism focuses on the axial anomaly and the special role played by the t]' 
in QCD [HIE], a detailed investigation of this state seems particularly desirable. The rj' requires 
disconnected contributions, and this poses a technical challenge. However, since the underlying 
physics is common to a broad class of vector-like gauge theories, there is no need to attack the 
problem in QCD. In this article the flavor-singlet state is studied in the generalized Schwinger 
model (QED in 2D, with Nf = 0, 1,2 massive flavors) [23j, which is much easier to simulate. 
In 2D a staggered field contains only two species. Hence, a square-root is required for Nf =1, 
while the Nf =2 continuum limit is supposed to be correct by definition. The point is that 
the conceptual issues match those of QCD. The rj in this model plays the same role as the t]' 
in QCD with three dynamical flavors, since its mass is predominantly due to the (global) axial 
anomaly. In the chiral limit of the Nf =1 theory it is known as the Schwinger particle. 



2 Simulation setup 

The goal is to perform t] spectroscopy with (rooted) staggered quarks in the massive Schwinger 
model {Nf =0, 1, 2) at several values of the coupling, such that the continuum limit (a— J-O) can 
be taken. In the Nf =1 case we wish to perform, in the second step, a chiral extrapolation to 
compare to the analytic prediction M^ = e^/7r at m = by Schwinger [23] . 

Because of the super-renormalizability of the Schwinger model [23], a convenient choice 
of scale is through the dimensionful coupling e in /3 = l/(ae)^. With this choice it is then 
straightforward to select the spatial extent Li = L, the temporal extent L2 = T and the coupling 
/3 such that eL is fixed (modulo cut-off effects). Moreover, due to the predictions of the eta 
mass in the chiral limit (M^ = i\^e^/7r for Nf =1,2, see [23]) and of the pion mass as a function 
of the quark mass (M^ = 2.008e^/^m^/^ for Nf =2, see [21]), one knows limm-s.o M^L beforehand, 
and one may choose the quark masses, at least for Nf =2, such that M-^^L assumes predefined 
values (again modulo cut-off effects). The parameters of the square lattices used in this article 
are shown in Tab.iH Most of them yield eL = eT = 17.89, but at one coupling (with five quark 
masses) dedicated finite- volume scaling studies are performed (see below). 

The covariant derivative in -Dgtag.m (and in the staggered 2-hop operators described below) 
uses 3 steps of APE smearing with a = 0.5 and back-projection to U(l). A great simplification 
is that reweighting techniques prove effective in 2D |25j. The plan is thus to generate substan- 
tial numbers of quenched lattices, and to include the determinant factor into the observable. 



GSM: p=12.8, L=64, mass=0.036, n_smear=3 
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Figure 1: Part of the Nj =0 Monte Carlo histories of the action s„ii (left) and of the 3- fold 
smeared field-theoretic topological charges ^'naijO'eeo (I'lght) on the 64^ lattices at /9 = 12.8. 



Using standard LAPACK routines, log(det(Dstag,m)) may be calculated via a Cholesky factor- 



stag,m-Dstag,m = R^ R by Summing logs of the diagonal elements of R. With R in hand, 
may be computed through two solve-for operations. Alternatively, one may use the LU 



IZatlOn 'LJstag,m-'-^sta,g,m 
-"-^stagjm 



decomposition with pivoting to compute the determinant and the inverse. 



3 Topological charge decorrelation 

Since the disconnected part of the rj correlator (for A^ = 1, 2) is sensitive to the global topological 
charge q of the gauge background U [261127]. it is crucial to achieve excellent ergodicity w.r.t. q. 
In the Schwinger model the electric flux of a gauge configuration is quantized; with toroidal 
boundary conditions [and ignoring a subset of configurations with measure zero] the quantity 



q[U] = — jd^xe^^Ff,^{x) = — J d^x F^ix) 



(2) 



takes on integer values. Field configurations with constant flux are known as instantons. They 
minimize the action in a given charge sector, but unlike in QCD they are completely delocalized. 
Such one-instanton configurations can be directly put on the lattice. Specifically [28j 



Ui{x) = exp(-27ria;2/[-^i-^2]) 
U2ix) = exp(+27rixi/Li ■ (Ja-^^iJ 



(3) 



is an implementation in pseudo- Coulomb gauge (where all links in the time direction are one, 
except those at X2 = L2). An anti-instanton follows by reversing the signs in the exponent. 

It is thus natural to realize a topology-changing update by multiplying a given configuration, 
link by link, with (|3]) or its conjugate (equal probability), the result being subject to a Metropolis 
accept-reject step. This way detailed balance is maintained by construction. With dynamical 
fermions, one would calculate det^-' (D[f/new])/det^^ (D[f/oid]) and include it into the decision. 
Upon combining this global update with a standard procedure [SUED], one has a simulation 



algorithm of a two-dimensionau U(l) gauge theory which may take large stepqj in configuration 
space. To further improve the symmetry of the overall charge distribution, one may perform, 
once in a while, a P-transformation of the gauge field, i.e. apply the transformation [31] 



f/,,(x,i)^f/^(-x-/i,t)^ (/i = l,..,rf-l), f/^(x,t)^f/^(-x,t) (/i = d) (4) 

which leaves the gauge action invariant but reverses the topological charge. 

In this paper, an A^ =0 gauge update consists of 4 over-relaxation sweeps [30] per Metropolis 
sweep [2S] (each Metropolis update applies 4 successive hits on a given link), and all of this is 
repeated L/{4a) times. In total, two adjacent configurations are separated as 

for n=l:L/4 

perform 2 over-relaxation sweeps 

perform 1 instanton/anti-instanton update 

perform 2 over-relaxation sweeps 

perform 1 Metropolis sweep with 4 hits per link 
end 
perform, with probability 0.5, a P-transformation 

and Fig.[T] shows a part of the pertinent Monte-Carlo history of the plaquette action and of 
the topological charge on the finest lattice (/3 = 12.8). The Wilson action per site is defined as 
Swi\{x) = l-ReUi2{x) = l-cos(6'i2) with Uuix) = Ui{x)U2{x + ei)Ul{x + e2)U^{x) = exp(i6'i2) 
denoting the plaquette. For the topological charge two field-theoretic definitional are used, 
Qnl = E sin(^!2V(2vr) G R ("naive") and gg), = E 0^^^/{2tt) G Z ("geometric"), where ^iJ 
denotes the plaquette angle after 3 smearing steps. The algorithm is seen to tunnel well. 

4 Staggered spectroscopy with all-to-all technology 

In weak coupling perturbation theory staggered fermions have been proven to be sensitive to the 
axial anomaly [32||33]. Still, it seems worth demonstrating that this sensitivity carries over, at 
the non-perturbative level, to the asymptotic states of the theory, and leads to a non- vanishing 
mass of the combined taste- and-fiavor singlet pseudoscalar meson in the chiral limit. 

The remnant staggered form of the continuum index theorem has been investigated in a 
classic paper [28]. A direct check of the mass excess of the rj {r]') in Nf =2 [Nf = 3) QCD 
with rooted staggered quarks and standard taste assignment has been attempted [3iH37] . but 
unfortunately in the disconnected contributions the signal dies quickly in the noise [381 - H3] . 
Encouraged by [2B1J2Z|, we now attack the same goal in the much simpler Schwinger model. 

With a single staggered field, we expect to find 4 pseudoscalar bosons. The lightest (750^5 ? 
to be dubbed 7r°) becomes massless (on an infinitely large lattice) in the limit m— J-O, the next 
two (75®G and 750^1^5, to be dubbed n"^) become massless up to cut-off effects, while the last 
one (75®!, to be dubbed rj) is supposed to be well separated and stay massive in the chiral 
limit. In standard terminology, the 7r° is taste-pseudoscalar, the vr^ are taste- (axial) vector (in 
2D there is no distinction), while the 77 is taste-scalar (or taste-singlet) - see [3J for details. 



^In fact, this idea may be used in the U(l) theory in 4D too, by just combining two such planar instantons 
(e.g. in the 12 and 34-planes, see [28] for details), albeit with the proviso that q then changes only by ±2 units. 

•^To avoid any correlation between adjacent configurations, on average a total drift by (xtop^)^ = io^)^ 
units must be realized, tantamount to a successful completion of approximately Xtop^ = ('7^) changes of q. 

''For q^^ normally a renormalization factor Z = l + 0{a^) is introduced, but with smearing this factor is so 
close to 1 that it seems permissible to drop it [in line with neglecting other 0{a^) efi^ects]. 
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GSM: p=7.2, L=48, mass=0.048, n_smear=3 
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• q_stag vs. q_geo 



GSM: |3=7.2, L=48, mass=0.048, n_smear=3 
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Figure 2: Left: Correlation between the topological charges q^l and qltlg, based on the T^q 
operator ^, with 3 smearings. Right: Overall topological charge distribution for Nf =0, and 
after reweighting to Nf =1, 2. Either plot refers to the intermediate mass run at f3 = 7.2. 



Staggered spectroscopy is performed by constructing dedicated correlators which project to 
a specific spinor-taste combination. In phenomenological apphcations it is common practice to 
use a source at a single lattice point as a "broad band emitter" which couples to all spinor-taste 
combinations, and to apply the projection only at the sink. In this work, we use an all-to-all 
propagator technique, and shall apply spinor-taste projection independently at the source and 
the sink. In explicit terms, for the 75(8^5 "Goldstone" state we use the operator Fss = (— l)^!"'"^^, 
which is a point-like operator. By contrast, the operator Fso — (75®!) is defined by [28j 
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r/.(a;,i/) 






— F2F1) where 



(5) 
(6) 



with ri^{x) = (— 1)^''<''^'', and is thus a 2-hop operator (it would be 4-hop in 4D). 

The operator (^ is supposed to be sensitive to the topological charge of the gauge back- 
ground U, and to test our implementation we routinely determine the fermionic chargq^ 



771 



qst>,g[U] = — tr(GF5o) 



(7) 



where G denotes the Green's function of the massive staggered operator. A scatter plot of (I7j) 
(y-axis) versus the gluonic charge qgeo[U] (x-axis) is shown in Fig.El We also checked that the 
operator F50F55 ~ leg) ,^5 is not sensitive to the topological charge. 

With these projections it is straightforward to construct the connected vr" and 7] correlators 

^ ^ G{x,y)T,,{y)G{y,x)T,,{x) (8) 



C^t) 
C,{t) 



L2 



x,y,X2—y2=t (mod T) 



^ E 



L2 



G{x, y')T5o{y', y)G{y, x')F5o(x', x) 



(9) 



x,y,X2—y2=t (mod T) 



For gg^g'g the statement in footnote 3 applies likewise (cf. Fig. [5] to see how close to 1 the Z- factor is). 
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Figure 3: aMesit) of Cj^it) (left) and Cr,{t) (right), at (3 = 7.2 for the intermediate quark mass. 



where 6{x',x)Tr,^{x) = T^r,{x',x) is used, and the primed positions are imphcitly summed over. 
For the combined taste-and-flavor singlet state there is also the disconnected contribution 



Dr,{t) = -^ E G{x,x')T,o{x',x)G{y,y')T,o{y',y) 

x,y,X2—y2=t (inodT) 

and for staggered fermions it must be combined with ([9]) in the form [2 

Fnit) ^ -^C,{t) - -LD.it) 



(10) 



N, 



m 



:iii 



to obtain the full 2-point function of the ri state. Here Nt = 2'^''^ denotes the number of "tastes" 
of a staggered field m d spacetime dimensions. Smee C,(() and f,(() fall off exponenfallfl 
at large t, with masses M^°°'^ and M^"" respectively (where only the latter one is physical), it 
follows that the ratio of the disconnected over the connected correlator takes the form 
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^ Z^ _ const e-^"* (12) 



N, 



f 



where zyVf^ = M* — M^°'^°, and the simplification applies for a<^t <^T. In other words, the 
prediction is that in 2D the ratio (IT2l) levels off at 2 for Nf =1, and at 1 for A^ =2. 



5 Telling Goldstone bosons from non-Goldstone bosons 

With the lattices of Tab.[T]in hand it is advisable to first check how well the different topological 
charges agree, how well the overall distribution is sampled, and whether the disconnected piece 
Drj(t) is indeed sensitive to the overall charge of the background U. 



^Strictly speaking, there is no transfer matrix argument ensuring this for the case of interest. For C^(t) one 
may be able to construct a transfer matrix in the partially quenched sense [44j . For F^ (t) , a transfer matrix exits 
only for Nf —2. For Nf ^l, there is no such argument (because of the rooted determinant), but our data are 
consistent with Dri{t) being the difference of two exponentials for Nf =0, 1, 2 alike. In other words, regardless of 
A^ the disconnected piece seems to have precisely the form needed to make pT|) a single exponential at t^oo. 
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Figure 4: Correlation between the disconnected Dn{t) and the squared topological charge q^ 
fort = T/4 (left) and t = T/2 (right). Data are from the intermediate quark mass run at 13 = 7.2. 

GSM: p=7.2, L=48, mass=0.016, n_smear=3 GSM: p=7.2, L=48, mass=0.080, n_smear=3 
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Figure 5: Ratio Rr^it) ( IT2j) at /3 = 7.2 for the smallest (left) and largest (right) quark mass. 



The correlation between any pair of the three topological charges considered is very good; 
for the larger two /3 even the Z-f actor inherent in q^^^^ and g^^^ is extremely close to 1, as can 
be inferred from Fig.[2j For Nf =0 the overall topological charge distribution is nicely sampled 
and close to a Gaussian. Under reweighting to A^ = 1 or A^ = 2 tiny asymmetries seem to get 
considerably enhanced. A feature relevant in what follows is that the dynamical distributions 
are narrower than the quenched one; with observables which are sensitive to topology it is 
useful to have the tails "oversampled" and to reduce their weight in the analysis. 

After these checks have been carried out successfully, it is straightforward to determine the 
pion mass aM.,^ for all couplings and quark masses. To this aim we consider the effective mass 



aM^it) = - log{[C{t-l) + ^Cit-iy-C{T/2y]/[Cit+l) + ^C{t+iy-CiT/2y]) (13) 





aM^ 


aMIrj 


/3 L/aNf 


mi 7712 "^3 "^4 "^5 


7771 7772 ^^5 '774 fn^ 


1.8 24 
1 
2 

3.2 32 
1 
2 

7.2 48 
1 
2 

12.8 64 
1 
2 


0.279(01) 0.368(01) 0.444(01) 0.513(01) 0.579(01) 
0.220(01) 0.318(01) 0.402(01) 0.478(01) 0.550(01) 
0.199(06) 0.292(02) 0.374(02) 0.454(01) 0.530(01) 
0.215(01) 0.281(01) 0.337(01) 0.389(01) 0.439(01) 
0.166(01) 0.240(01) 0.304(01) 0.362(01) 0.417(01) 
0.154(07) 0.217(04) 0.283(01) 0.342(01) 0.400(01) 
0.142(01) 0.186(01) 0.224(01) 0.260(01) 0.294(01) 
0.111(02) 0.159(01) 0.203(01) 0.242(01) 0.279(01) 
0.091(16) 0.137(06) 0.188(02) 0.230(01) 0.269(01) 
0.107(01) 0.140(01) 0.169(01) 0.196(01) 0.221(01) 
0.087(02) 0.123(01) 0.153(01) 0.183(01) 0.210(01) 
0.086(07) 0.114(02) 0.143(02) 0.175(01) 0.201(02) 


0.39(02) 0.36(07) 0.33(03) 0.31(04) 0.30(05) 
0.83(42) 1.05(60) 0.63(08) 0.57(14) 0.53(21) 

0.28(02) 0.22(01) 0.21(01) 0.21(01) 0.18(01) 
1.41(97) 0.48(02) 0.39(06) 0.36(02) 0.40(06) 

0.17(02) 0.14(01) 0.12(01) 0.11(01) 0.12(03) 
0.45(07) 0.29(02) 0.28(01) 0.24(02) 0.22(02) 

0.14(02) 0.12(01) 0.10(02) 0.09(01) 0.09(02) 
0.32(05) 0.24(04) 0.23(02) 0.19(01) 0.17(02) 


3.2 24 
1 
2 

3.2 40 
1 
2 


0.231(01) 0.285(01) 0.338(01) 0.390(01) 0.439(01) 
0.182(01) 0.246(01) 0.306(01) 0.363(01) 0.416(01) 
0.173(02) 0.229(01) 0.288(01) 0.346(01) 0.400(01) 
0.214(01) 0.280(01) 0.337(01) 0.389(01) 0.439(01) 
0.158(02) 0.241(01) 0.302(01) 0.363(01) 0.417(01) 
0.138(04) 0.220(04) 0.276(05) 0.347(02) 0.402(01) 


0.26(01) 0.23(01) 0.21(01) 0.19(01) 0.18(01) 
2.12(1.56) 0.64(16) 0.47(02) 0.37(04) 0.33(05) 

0.28(13) 0.26(03) 0.21(02) 0.18(03) 0.21(04) 
2.75(2.72) 1.23(83) 0.41(25) 0.38(12) 0.36(04) 



Table 2: Measured qMt^ (Goldstone boson) and aAM^ (excess of the taste-and-Havor singlet 
pseudoscalar state) for all /3, L/a, Nf , and m. The quark masses are given in Tab. 1. 



with C = Ctt, which is designed to compensate for the influence of the periodic boundary 
conditions. A typical example is shown in Fig.|3l A nice plateau is observed; albeit with tiny 
oscillationqj which grow towards the center of the box. Under reweighting to Nf =1 or Nf =2 
this effect gets enhanced. For the connected piece of the eta (which is an unphysical state) 
similar results are obtained. Throughout, the difference between these two masses is tiny. 

Turning to the disconnected piece /^^(t), we first consider its correlation with the topological 
charge of the background. Example results are presented in Fig. HI For these t the sub-ensemble 
average {D^[t))q seems to be a linear function of g^ (which, from a glimpse at flTOj) and the 
definition of Q'stag[^] is plausible). The issue most relevant is whether reasonable results for the 
disconnected-over-connected ratio ( TT2l) are obtained. Typical results are presented in Fig.O 
We obtain a rather clear signal up to about t = T/4 and find a qualitatively different behavior 
for the three Nf shown. In the quenched case the pattern is consistent with a linear rise (with 
a slope which clearly depends on the quark mass) . After reweighting to Nf =1 or Nf =2 the 
behavior is consistent with the prediction flT2l) . A typical problem with the disconnected piece 
is that at large t the data may go astray without the statistical error giving a hint of this (we 
tried several jackknife blocksizes). Still, we can extract the mass gap a/SM^ from a fit to (1121) at 
intermediate t, i.e. before the noise prevails. The results such obtained are collected in Tab. [21 
along with the aM^, values determined previously. 

The lower part of Tab. [2] contains the results of a finite-volume scaling study at /3 = 3.2 
(since finite volume effects relate to infrared physics the restriction to one /3 is permissible). 
It seems that in the extra-small volume the pion mass is affected for the lightest two quark 
mass values. In the standard volume only the lightest quark mass data suffer from small finite 
volume artefacts. In case of the mass gap aAM^ no finite volume effects are found ioi Nf =1, 



^This may signal the presence of a parity partner [3] on odd timeslices. In our analysis we use M^(T/2 — 1) 
which, in turn, is based on the correlator C^(i) on the even timeslices t = r/2 — 2,r/2. 
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while for Nf =2 the quahty of the data is less convincing, which likely indicates the limitations 
of the reweighting method (fortunately, we are mostly interested in aAM^ for Nf =1). 

The next step is a continuum extrapolation of the meson masses obtained. We do this both 
for Mtt/c in the Nf =2 theory (to test the prediction by Smilga [21], see above) and for the 
main object of interest AM^/e in the Nf =1 theory. As mentioned in Sec. 2 the simulations 
were carried out at fixed physical quark masses (i.e. m/e), such that the data of Tab. |2] can be 
continuum extrapolated without interpolation in m. These extrapolations are shown in Fig.|6l 
It seems that all four lattice spacings are in the Symanzik scaling regime; we obtain acceptable 
linear fits with two degrees of freedom. 

With these continuum results for M^^/e (for Nf =2) and AM^/e (for Nf = 1) in hand, it is 
interesting to consider their quark mass dependence. In the left panel of Fig. [7] we plot {M.,^/e)^ 
as a function of (m/e)^. Here, the lightest pion mass has been adjusted by the finite volume 
correction factor 0.138/0.154 = 0.896 (taken from Tab. [2]). Using a 1-parameter ansatz fits the 
data with x^/d.o.f = 1.13 and a slope parameter of 8.221(63). This is in reasonable agreement 
with Smilga's prediction that this parameter should be 2.008^ = 8.096 ^^. 

Last but not least let us consider the quark mass dependence of the mass gap AM^/e in the 
Nf =1 theory, as shown in the right panel of Fig. [71 Being unaware of an analytic prediction, 
it is not entirely clear which powers of AM^/e and m/e one should choose to display the 



data, and we opt for staying without additional powers. Fortunately, in this representation 
the quark mass dependence seems mild, and both a linear and a quadratic fit with one degree 
of freedom describe the data convincingly. Taking half of the spread as the systematic error, 
these fits predict AM^/e = 0.591 (25) (37) in the chiral limilo, which is in perfect agreement 
with Schwinger's analytic result M^/e = AM^/e = l/\/7r = 0.56419 [23] . 

6 Summary 

Triggered by the criticism of |1H0] the validity of the "rooting trick" in studies with staggered 
quarks has been the subject of an intense debate in the lattice community. Several review talks 
at major conferences [71-[T0] presented a wealth of numerical and analytical evidence in favor of 
the procedure, but so far the "experimentum crucis", i.e. a direct test of r^'-phenomenology in 
QCD with rooted staggered quarks remained elusive, due to noise issues [5iHS7] . 

This paper is based on the observation that there is no strict need to investigate the topic 
in QCD, since the conceptual issue is one-to-one matched in the massive Schwinger model with 
1 fiavor, which is much simpler to simulate. We demonstrated that in this case it is possible to 
obtain conclusive results for the disconnected-over-connected ratio f lT2|) . As shown in Fig. [5] the 
prediction flT2|) with Nf =1 and Nt = 2 is beautifully confirmed. In consequence the mass gap 
AM^i or the physical mass M^ in the 1-fiavor theory can be determined with sufficient precision, 
so that a combined continuum and chiral extrapolation is possible. The result is in perfect 
agreement with the analytical prediction M^ = e/y^ by Schwinger |23j . 

It seems this is the first ad oculos demonstration that the staggered setup - with the rooting 
trick in the functional measure - treats the contribution of the axial anomaly to the particle 
spectrum correctly. With this result, and in view of [51HTF] . one may predict that the outcome 
of a similar study in QCD will be identical, once sufficient CPU power is available. 
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